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ABSTRACT 

We  develop  and  describe  a  stable  gain  scheduling  controller 
for  a  gas  turbine  engine  that  drives  a  variable  pitch  propeller.  A 
stability  proof  is  developed  for  gain  scheduled  closed-loop  sys¬ 
tem  using  global  linearization  and  linear  matrix  inequality  (LMI) 
techniques.  Using  convex  optimization  tools,  a  single  quadratic 
Lyapunov  function  is  computed  for  multiple  linearizations  near 
equilibrium  and  non-equilibrium  points  of  the  nonlinear  closed- 
loop  system.  This  approach  guarantees  stability  of  the  closed- 
loop  gas  turbine  engine  system.  Simulation  results  show  the  de¬ 
veloped  gain  scheduling  controller  is  capable  of  regulating  a  tur¬ 
boshaft  engine  for  large  thrust  commands  in  a  stable  fashion  with 
proper  tracking  performance. 

NOMENCLATURE 

N\  Low  Spool  Speed  (Non-dimensional) 


N2  High  Spool  Speed  (Non-dimensional) 
u\  Fuel  Flow  Control  Input  (Non-dimensional) 

M2  Propeller  Pitch  Angle  Control  Input  (Deg) 

T  Thrust  (N) 

TSFC  Thrust  Specific  Fuel  Consumption  (kg/s/N) 
a  Scheduling  Parameter 
A  Eigenvalue 


1  Introduction 

Stability  and  control  of  gas  turbine  engines  have  been  of  in¬ 
terest  to  researchers  and  engineers  from  a  variety  of  perspectives. 
An  introduction  to  the  analysis  and  design  of  engine  control  sys¬ 
tems  can  be  found  in  [1],  The  basics  of  controlling  a  gas  turbine 
engine  while  satisfying  numerous  constraints  has  been  reviewed 
in  [2],  The  design  of  engine  control  and  monitoring  systems  with 
a  dual  interest  in  both  turbofan  and  turboshaft  engines  has  been 
covered  in  [3].  An  application  of  robust  stability  analysis  tools 
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for  uncertain  turbine  engine  systems  is  presented  in  [4].  An  ap¬ 
plication  of  the  Linear-Quadratic-Gaussian  with  Loop-Transfer- 
Recovery  methodology  to  design  a  control  system  for  the  F-100 
turbofan  engine  is  presented  in  [5],  and  for  a  simplified  turbofan 
engine  model  is  considered  in  [6].  A  unified  robust  multivari¬ 
able  approach  to  propulsion  control  design  has  been  developed 
in  [7].  The  development  of  other  control  techniques,  such  as 
sliding  mode,  for  gas  turbine  engine  application  can  be  found 
in  [8]. 

To  facilitate  the  stability  analysis  of  nonlinear  systems,  such 
as  gas  turbine  engines,  an  efficient  technique  is  to  approximate 
them  by  a  linear  time-varying  (LTV)  system.  This  concept, 
which  is  known  as  global  linearization,  can  be  found  in  [9, 10]. 
More  recent  work  on  global  linearization  and  the  use  of  Linear 
Matrix  Inequalities  (LMIs)  for  the  analysis  of  dynamical  sys¬ 
tems  can  be  found  in  [11].  Some  Soviet  literatures  on  the  ab¬ 
solute  stability  problem,  like  Lur’e  and  Postnikov  [12, 13]  and 
Popov  [14-16],  also  implicitly  use  the  idea  of  global  lineariza¬ 
tion.  Recent  literatures  that  demonstrate  the  practical  power  of 
global  linearization  technique  include  [17,  18],  and  [19].  [19] 
uses  this  idea  along  with  the  notion  of  incremental  stability. 

One  of  the  control  design  approaches,  which  perhaps  is  one 
of  the  most  popular  nonlinear  control  design  approaches  and 
has  been  widely  and  successfully  applied  in  fields  ranging  from 
aerospace  to  process  control  [20,  21],  is  gain  scheduling.  Gas 
turbine  engines  are  no  exception,  and  research  on  gain  schedul¬ 
ing  control  of  gas  turbine  engines  is  presented  in  [22-29].  A 
simplified  scheme  for  scheduling  multivariable  controllers  for 
robust  performance  over  a  wide  range  of  turbofan  engine  op¬ 
erating  points  is  presented  in  [23].  In  a  recent  work  presented 
in  [27],  results  on  polynomial  fixed-order  controller  design  are 
extended  to  SISO  gain- scheduling  with  guaranteed  stability  and 
Hoo  performance  for  a  turbofan  engine,  over  the  whole  schedul¬ 
ing  parameter  range.  In  this  work,  the  engine  Linear  Parameter 
Varying  (LPV)  representation  depends  on  an  exogenous  variable 
parameter  which  is  the  combustion  chamber  pressure. 

In  the  previous  works  [30,  31]  the  authors  discussed  con¬ 
trollers  for  single  spool  and  twin  spool  turboshaft  systems.  Those 
controllers  are  designed  for  small  transients,  and  small  throttle 
commands.  In  this  paper  we  develop  an  output  dependent  gain 
scheduled  control  structure  for  a  MIMO  linear  parameter  depen¬ 
dent  model  of  a  JetCat  SPT5  turboshaft  engine  [32]  using  the 
method  presented  in  [20,33-35].  This  controller  is  designed  to  be 
used  for  the  entire  flight  envelope  of  the  twin  spool  turboshaft  en¬ 
gine  with  stability  guarantees.  The  scheduling  variable  in  our  de¬ 
sign  process  is  an  endogenous  parameter,  which  is  a  function  of 
the  gas  turbine  engine  spool  speeds.  This  endogenous  scheduling 
variable  captures  the  plant  nonlinearities,  as  explained  in  [33,34], 
since  the  spool  speeds  are  the  main  states  of  the  turboshaft  engine 
state-space  model,  and  also  the  outputs  of  the  system.  The  stabil¬ 
ity  analysis  for  the  closed-loop  system  presented,  and  the  essen¬ 
tial  part  of  the  stability  analysis  is  to  find  a  common  quadratic 


Lyapunov  function  for  multiple  linearizations  near  equilibrium 
and  non-equilibrium  points,  which  are  distributed  over  the  entire 
operational  envelope  of  the  plant. 

The  paper  is  organized  as  follows:  First,  a  linear  parameter 
dependent  representation  of  the  plant  is  presented.  Concepts  for 
output  dependent  gain  scheduled  control  of  this  model  are  then 
developed.  Third,  the  stability  analysis  of  the  closed-loop  sys¬ 
tem  is  presented.  Finally,  simulation  results  for  gain  scheduling 
control  of  a  MIMO  physics-based  model  of  a  JetCat  SPT5  tur¬ 
boshaft  engine  are  presented.  Simulation  results  show  the  suc¬ 
cessful  application  of  the  proposed  controller  for  the  entire  flight 
envelope  of  the  turboshaft  engine  with  guaranteed  stability  and 
proper  tracking  performance. 


2  Gain  Scheduling  Control  Design 

Consider  the  nonlinear  dynamical  system 

xP  =  fP(xP,u), 

(1) 

y  =  gP(xP,u ), 

where  xp  G  is  the  state  vector,  u  G  is  the  control  input 
vector,  y  G  is  the  output  vector,  fp(.)  is  an  ^-dimensional  dif¬ 
ferentiable  nonlinear  vector  function  which  represents  the  plant 
dynamics,  and  gp(.)  is  an  m- dimensional  differentiable  nonlinear 
vector  function  which  generates  the  plant  outputs.  We  intend  to 
design  a  feedback  control  such  that  y  properly  tracks  a  reference 
signal  r  as  time  goes  to  infinity,  where  r  E  Dr  C  9tm,  and  Dr  is  a 
compact  set. 

Assume  that  for  each  r  G  Dr ,  there  is  a  unique  pair  (x%,ue) 
that  depends  continuously  on  r  and  satisfies  the  equations 

0  =  fp{4,ue), 

(2) 

r  =  gp(xpe,ue ), 

xp  is  the  desired  equilibrium  point  and  ue  is  the  steady-state  con¬ 
trol  that  is  needed  to  maintain  equilibrium  at  xp.  It  is  often  useful 
to  parameterize  the  family  of  system  equilibria  as  follows: 

Definition  1.  The  functions  xp  (a)  ,ue  (a),  and  re(a)  define  an 
equilibrium  family  for  the  plant  (1)  on  the  set  Q  if 

fp(xp(a)Ma))  =  0, 

gp(x%(a),ue(a))  =  re(a),  ae£2. 

Let  O  C  be  the  region  of  interest  for  all  possible  sys¬ 
tem  state  and  control  vector  (. xp,u )  during  the  system  operation, 
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and  denote  xpei  and  ueu  i  E  I  =  1,2,  ...,g,  as  a  set  of  constant  op¬ 
erating  points  located  at  some  representative  and  properly  sepa¬ 
rated  points  inside  O.  Introduce  a  set  of  q  regions  O*-,  i  E  /  cen¬ 
tered  at  the  chosen  operating  points  (xpei^uei),  and  denote  their 
interiors  as  O/o,  such  that  Oyo  Pi  O&o  =  0  for  all  j  /  k,  and 
U-=1  Q  =  O.  The  linearization  of  the  plant  at  each  equilibrium 
point  is 


xp  =  Ap  (xp  -  Xpei )  +  Bp  (u  -  uei ) , 

V  =  cf  (xp  -  x£.)  +  Df  (u  -  Uei)  +  yei, 


where  the  matrices  are  obtained  as  follows 


Af  = 

BPi  = 

Cf  = 

Df  = 


df_ 

dxP 

dfP_ 

du 

dgP_ 

dxP 

dgP 


0 ei,uei)  5 

\/(xp,u)  e  Oi , 

(■ 5 

V(xp,u)  e  Oi, 

1(4^)’ 

\/(xp,u )  e  Oi, 

{^ei^ei)  ’ 

"i(xp  ,ii)  e  o,. 

(4) 


(5) 


model.  Then,  in  between  operating  points,  the  controller  gains 
are  linearly  interpolated  such  that  for  all  frozen  values  of  the 
parameters,  the  closed-loop  system  has  satisfactory  properties, 
such  as  nominal  stability  and  robust  performance.  To  guaran¬ 
tee  that  the  closed-loop  system  retains  the  dynamic  properties 
of  the  frozen-parameter  designs,  the  scheduling  variables  should 
vary  slowly  with  respect  to  the  system  dynamics  [33].  Figure  1, 
shows  schematically  how  the  output  dependent  gain  scheduled 
controller  works. 


Scheduling 


FIGURE  1 .  Output  dependent  gain  scheduling  controller  diagram 


Note  that  (xp,u)  belongs  to  only  one  O*  at  each  time.  Corre¬ 
sponding  to  each  linearization  at  ith  equilibrium  point,  there  ex¬ 
ist  an  at  E  £2,  which  is  a  function  of  equilibrium  values  of  the 
system  outputs,  i.e.  ye\. 

The  family  of  plant  linear  models  (4)  can  be  written  as 

8xp  =Ap(a)8xp+Bp(a)8u, 

(6) 

Sy  =  Cp(a)Sxp  +  Dp(a)Su ,  Va  e  £2, 

where 

8xp  =  xp  —  xp(a) 

8y  =  y-ye{a),  (7) 

8  u  =  u  —  ue(a). 

Ap(a),  Bp(a),  Cp(a),  and  Dp(a)  are  the  parameterized  plant 
linearization  family  matrices  and  xp (a),  ue(a),  and  ye{oc)  are 
the  parameterized  steady-state  variables  for  the  states,  inputs  and 
outputs  of  the  plant,  which  form  the  equilibrium  manifold  of 
plant  (1).  The  subscript  ’V’  stands  for  ’’steady- state”  through¬ 
out  this  paper. 

Based  on  the  results  from  [20,  33-35],  an  output  depen¬ 
dent  gain  scheduled  controller  for  plant  (6)  is  designed  as  fol¬ 
lows.  First,  a  set  of  parameter  values  a;  is  selected,  which 
represent  the  range  of  the  plant’s  dynamics,  and  a  linear  time- 
invariant  controller  is  designed  for  each  corresponding  linear 


The  parameter  a  is  called  the  scheduling  variable  and  should 
be  measurable  in  real  time,  a  can  be  a  function  of  endogenous 
variables  (i.e.,  depending  on  the  plant  states)  and/or  exogenous 
variables  (i.e.,  independent  of  the  plant  states).  In  LPV  systems, 
this  parameter  is  an  exogenous  parameter  [36].  Some  of  the  ex¬ 
amples  of  exogenous  parameter  selection  in  LPV  control  of  tur¬ 
bine  engines  are  presented  in  [25-27].  In  [25],  the  scheduling 
parameter  is  defined  as  a  function  of  the  exogenous  signals  de¬ 
scribing  the  surroundings,  like  altitude,  intake  Mach  number,  and 
a  health  parameter  describing  the  state  of  the  compressor.  In  [26], 
the  scheduling  parameter  is  defined  as  a  function  of  lagged  mea¬ 
surement  of  engine  thrust  and  altitude,  which  are  exogenous  vari¬ 
ables.  In  [27],  the  scheduling  parameter  is  defined  to  be  the  com¬ 
bustion  chamber  pressure,  which  is  an  exogenous  variable.  In 
gain-scheduling,  this  parameter  is  a  function  of  the  output  and 
hence  it  is  an  endogenous  parameter  [36].  Some  of  the  examples 
of  endogenous  parameter  selection  for  gain- scheduled  control  of 
turbine  engines  can  be  found  in  [22,28,29].  In  [22],  the  schedul¬ 
ing  parameter  is  defined  to  be  the  engine  low  pressure  spool 
speed,  which  is  one  of  the  outputs  of  the  system.  In  [28,29],  the 
scheduling  parameter  is  defined  to  be  the  engine  high  pressure 
spool  speed.  In  the  turboshaft  engine  control  example  described 
later  in  this  paper,  a  is  defined  to  be  the  Euclidean  norm  of  the 
engine  spool  speeds,  which  can  be  measured  in  real-time.  Since 
the  spool  speeds  are  the  only  plant  states  in  the  model  and  also 
due  to  the  fact  that  we  need  the  plant  nonlinearities  to  be  cap¬ 
tured  by  the  output  vector,  as  explained  in  [33,  34],  we  defined 
the  scheduling  parameter  to  be  the  function  of  both  spool  speeds. 
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On  the  other  hand,  for  a  simpler  interpolation  process,  we  wanted 
the  scheduling  parameter  to  be  a  scalar,  so  we  used  the  Euclidean 
norm  of  the  output  vector. 

The  design  of  a  linearization  gain  scheduled  controller  re¬ 
quires  designing  a  linear  controller  family  corresponding  to  the 
plant  linearization  family  (6).  Let  the  parameterized  linear  con¬ 
troller  family  be 

Sxc  =  Ac  (a)Sxc  +  Bc  (a)[<5y  —  <5r], 

8u  =  Cc(a)Sxc  +  Dc(a)[Sy  —  <5r],  Va 


where 


Ac(a)  = 
Bc(a)  = 
Cc(oc)  = 
Dc(a)  = 


d_r_ 

dxc 

dfc 


dxc 

dgc 

dy 


l(j£(a),ye(a),re(a))> 
I  (x?e{a),ye(a),re{a))i 
I  (x%(a),ye(a),re(a))> 
\{xce{cc),ye(a),re{a)), 


Va  e  £1 


where 


So  the  controller  family  has  the  form 


(13) 


Sxc  =xc-xce(oc), 

Sr  =  r  —  re(a),  Va  G  Q. 


(9) 


xce(a)  and  re(a)  are  the  parameterized  steady- state  variables  for 
the  controller  states  and  reference  signals.  A  standard  realization 
of  the  parameterized  controller  can  be  written  with  the  reference 
signal  explicitly  displayed  as 


Ac(a)  Bc(a)  -Bc(a) 
Cc(a)  Dc(a)  - Dc(a ) 


8xc 

Sy 

- 

Sr 

VaeQ.  (10) 


We  have  to  obtain,  based  on  the  linear  controller  family  (10), 
a  controller  that  has  the  general  form 


xc  =  Ac(a)[xc  —  xce(a)\  -\- Bc (a)[y  —  r\, 

u  =  Cc(a)[xc  —  xce(a)]  +Dc(a)[y  —  r]  +  ue(a ),  Va  e  Q. 

(14) 

Note  that  re(a)  =ye  (a),  as  a  result  8y  —  8r  =  y  —  r.  The  schedul¬ 
ing  parameter  a  is  treated  as  a  parameter  throughout  the  de¬ 
sign  process,  and  then  it  becomes  a  time-varying  input  signal 
to  the  gain- scheduled  controller  implementation  through  the  de¬ 
pendence  a(t)  =  p(y(t)).  The  parameter  a  is  an  endogenous 
variable,  since  it  is  a  function  of  the  plant  outputs.  Replacing  a 
with  p(y),  the  gain  scheduled  controller  becomes 


Xc  =  Ac(p(y))[xc -X?e{p(y))}  +Bc(p(y))\y  -  r], 

u  =  Cc (p(y))[xc  - x?e(p(y))}  +Dc(p(y))\y-r\+ue(p(y)). 

(15) 

Linearization  of  (15)  about  an  equilibrium  specified  by  a  yields 


ic  =  /c(xc,j,r), 
u  =  gc{xc,y,r), 


(ID 


with  the  input  and  output  signals  corresponding  to  the  nonlin¬ 
ear  plant  (1).  fc(.)  is  an  m- dimensional  differentiable  nonlinear 
vector  function  which  represents  the  controller  dynamics,  and 
gc(.)  is  an  m-dimensional  differentiable  nonlinear  vector  func¬ 
tion  which  generates  the  controller  outputs. 

The  objective  in  linearization  scheduling  is  that  the  equilib¬ 
rium  family  of  the  controller  (11)  match  the  plant  equilibrium 
family,  so  that  the  closed-loop  system  maintains  suitable  trim 
values,  and  the  linearization  family  of  the  controller  obtained 
from  linearizing  (1 1)  is  the  same  as  the  designed  family  of  linear 
controllers  shown  in  (8)  [20].  Lor  the  equilibrium  conditions  of 
plant  (1)  and  controller  (1 1)  to  match,  there  must  exist  a  function 
x^(a)  such  that 


§xc  =  Ac (a)Sxc  +  Bc (a)[y -  r\ 

-[A^a)^^}x[^(ye(a))8y], 

Su  =  Cc{a)5xc  +  Dc(a)[y-r] 

,due{a) 


(16) 


da 


a  off]  x  [|fc(c ow. 


Comparing  (16)  with  (10),  we  see  there  are  additional  terms,  and 
we  refer  to  them  as  hidden  coupling  terms  following  the  notation 
of  [20].  In  order  to  get  rid  of  these  hidden  coupling  terms,  the 
following  condition  must  be  satisfied 


[At(a)^|^lx[^(j,(ci))Svl=0, 

|^>*C'(O)^]x[|(v«(c<))H  =  0. 


(17) 


It  is  not  always  easy  to  come  up  with  solutions  to  satisfy  condi¬ 
tion  (17).  In  order  to  make  the  design  process  easier,  we  control 
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0  =  /c«(a),)v(a),re(a)), 

ue(a)  =  gc{xce(a),ye(a),re(a)),  Va  e  Q., 


(12) 


the  system  via  filtered  inputs,  rather  than  the  input  themselves, 
so  there  is  no  need  for  equilibrium  control  value  other  than  zero 
(i.e.  *£(a)  =  0,  ve(a)  =  0,  Va,  where  ve(a)  is  the  parameterized 
steady-state  variables  for  the  new  inputs). 

The  plant  (1)  with  the  filtered  inputs,  becomes 


r  p 


F(xP,u ) 
-r\cu 


0 

ricxl 


(18) 


y  =  gp(Xp,u), 

The  controller  has  the  general  form 

x<  =  r(xc,y,r), 
v  =  gc(xc,y,r ), 


(19) 


with  the  input  and  output  signals  corresponding  to  those  of  the 
nonlinear  plant  (18).  Now,  combining  (18)  and  (19)  leads  to 


xp~ 

1 

0 

u 

= 

-Bcu 

+ 

t)c  X  / 

xc 

_r(xc,gp(xp,u),r)_ 

0 

*  fM 

v  =  gc(xc,gp(xp,u),r), 

S - V - ' 

gM 

and  the  closed-loop  nonlinear  system  is 

x  =  f(x,r)+Bg(x,r), 
=  F(x,r), 


(20) 


(21) 


where  x  G  Dx  c  9i''+2"',  and  r  £  D,  C  91'”.  The  augmented  linear 
family  of  systems  for  the  augmented  plant  (18)  becomes 


8xp 

\Ap(a)  Bp(a)  1 

8xp 

8u 

y  s 

O 

1 

C5 

X 

V 

8u 

5xaug 

- V* - 

^aug(o^) 

^JCaug 

8y=[Cp(a),Dp(a)\ 

s - V - ' 

Qtug(°0 


8xp 

8u 


(22) 


and  the  controller  has  the  similar  structure  as  (8) 

8xc  =  Acv(a)Sxc  +  Bcv(a)[Sy  —  <5r], 

<5v  =  Cy(a)8xc  +Z)J(a)[5y  —  <5r],  Va  e  £2, 


(23) 


where 

Sv  =  v  —  ve(a),  Va  e  Q. 

Now,  since  xce(a)  =  0,  ve(a)  =  0,  Va,  the  controller  is 

xc  =  Acv(a)xc  +  Bcv(a)[Sy  -  <5r], 
v  =  Cy(oc)xc  -\-D^{a)[8y  —  Sr],  Va  e  Q, 

rewriting  controller  (25)  with  a(t)  =  p(y(t )),  we  obtain 

xc  =  Acv  0 p(y )  )xe  +  (p(y) )  b  -  r] , 
v  =  Ccv(p{y))xc +Dcv{p{y))[y-r\. 


(24) 


(25) 


(26) 


Linearization  of  (26)  about  an  equilibrium  specified  by  a  gives 
(25),  so  there  are  no  hidden  coupling  terms  similar  to  the  ones  we 
saw  in  (16),  and  the  condition  (17)  is  satisfied.  One  of  the  options 
for  control  design  is  to  set  the  controller  matrices  as  follows 


Acv{a)  =ACV  =  —£CI,  Bcv(a)  =  BC  =  I, 
Cy(a)  =  Kj(a),  Dcv(a)  =  Kp(a), 


(27) 


which  is  a  kind  of  proportional-plus-integral  (PI)  control,  where 
Ki(a)  is  the  integral  control  gain  matrix,  and  Kp(a)  is  the  pro¬ 
portional  control  gain  matrix.  Hence  the  controller  for  the  aug¬ 
mented  plant  linearization  family  (22)  has  the  final  form 


-£CI  I  -I 

Ki{a)  Kp(a)  -Kp(a) 


V 

Sy 

- 

Sr 

Va  g  Q.  (28) 


The  linearized  closed-loop  system  (22)  with  controller  (25)  be¬ 
comes 


~  8xp~ 

Ap(a) 

Bp(a) 

0 

8ii 

= 

ricD^,(a)Cp(a) 

-ricl+Dcv(a)Dp(a) 

VcCy(a) 

xc 

_  Bcv(a)Cp(a) 

Bcv(a)Dp(a ) 

Acv(a) 

Sx 


8xp~ 

0 

8u 

+ 

~BcDcv(a) 

xc 

_  ~Bcv{a)  _ 

Ad(a) 

8r ,  Va  G  Q.. 


8x 


Bd  («) 


(29) 

For  the  case  where  we  have  plant  states  as  the  outputs  8y  =  8xp, 
(i.e.  Cp(a)  =  I,Dp(a)  =  0)  the  linearized  closed-loop  system 
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(22)  with  controller  (28)  becomes 


~8xp  1  [  Ap(a)  Bp (a)  0  1  [SxP~ 

Su  =  ricKp(a)  -r\cI  rfcKi(a)  8u 

_xc  \  [  1  0  -£CI  \  [  xc  _ 


Bd(oc) 


3  Stability  Analysis 

In  this  section  we  show  the  stability  of  the  closed-loop  non¬ 
linear  system  by  using  ’’global  linearization”  technique.  The  sta¬ 
bility  is  due  to  the  existence  of  a  single  quadratic  Lyapunov  func¬ 
tion  for  all  a  G  £2,  by  computing  a  single  Lyapunov  matrix  P 
using  convex  optimization  tools. 

Assumption  1 .  The  matrices  Aci  and  Bci  are  hounded 

||Ac/(0II<*a,  IIMOII<*fl,  Vr>0,  (31) 

where  and  ks  are  constants. 

To  analyze  the  stability  of  the  nonlinear  closed-loop  system,  we 
use  a  technique  known  as  ’’global  linearization”  developed  in  [9- 
11]. 

Theorem  1.  Consider  the  closed-loop  system  (21),  and  as¬ 
sume  there  is  a  family  of  equilibrium  points  (xe ,  re )  such  that 
F(xe,re )  =  0.  Define  A1^  =  £  S,  \/x  G  Dx,  where  S  is  the  set 

of  linearizations  of  system  (21) 

S:={A$,VxeDx}.  (32) 

Assume  there  exists  a  symmetric  positive  definite  matrix  P,  such 
that 

PA"l  +AnJlTP  <  0,  \/Anc{  G  S,  (33) 


where  £  S,  for  all  i  G  {1,2,...,L}.  Note  that  A*\. 

I X — Xf 

can  be  obtained  by  linearizing  the  nonlinear  system  (21)  at  non¬ 
equilibrium  points  (transient  condition),  and  also  at  equilibrium 
points  (steady  state  condition),  which  in  this  paper,  are  repre¬ 
sented  by  Aci(ai).  Then  using  convex  optimization  tools  [37, 38], 
we  compute  a  common  symmetric  positive  definite  matrix  P,  such 
that 

PA"lh+AfP<  0,  Vi G  {1,2, (35) 

With  assumption  1  satisfied,  we  can  verify  Acfia)  G  S,  for  all 
a  G  Q,  then  system  (29)  is  also  asymptotically  stable. 


4  Turboshaft  Engine  Example 

We  apply  the  proposed  output  dependent  gain  scheduling 
controller  to  a  physics-based  model  of  the  JetCat  SPT5  tur¬ 
boshaft  engine  driving  a  variable  pitch  propeller  developed  in 
[39,40].  Note  that  some  of  the  plant  states  and  inputs  have  been 
non-dimentionalized  by  their  design  values;  fuel  flow  input,  u\, 
is  divided  by  0.0035323  (kg/s),  core  spool  speed,  Afc,  which  is 
the  first  plant  state  (x\),  and  is  divided  by  170000  RPM,  and 
fan  spool  speed,  N\,  which  is  the  second  plant  state  (x^),  and  is 
divided  by  7000  RPM.  The  complete  setup  of  the  JetCat  SPT5 
turboshaft  engine  with  a  variable  pitch  propeller  on  the  test  stand 
is  shown  in  Figure  2. 

4.1  Equilibrium  Manifold 

For  a  standard  day  at  sea  level  condition  we  chose  five  prop¬ 
erly  separated  equilibrium  points  on  the  plant  equilibrium  man¬ 
ifold  for  linearizing  the  plant  model  at  those  points.  The  lin¬ 
earization  matrices  for  these  five  equilibrium  points  and  steady 
state  values  of  the  engine  variables,  scheduling  parameter  and 
control  parameters  are  given  as  follows: 

Equilibrium  Point  1  (Full  Thrust): 

ulel  =  1.0,  u2ei  =  16  (deg),  xui  =  1.0,  x2ei  = 

0.9524,  Tei  =  255.8685  (N),  a,  =  1.3810,  and 


then  the  system  (21)  is  asymptotically  stable. 

Remark  1 .  In  practice  we  can  not  obtain  S,  instead,  we  can 
linearize  system  (21  )for  a  large  number  of  states  Xi,  i  =  1, . . .  ,L, 
which  we  claim  is  sufficient  to  cover  the  set  of  actual  operating 
conditions,  to  show  the  stability  of  the  closed-loop  system.  Define 
S  as  a  matrix  polytope  described  by  its  vertices 


S.=  C  o{A£,„.,A*}, 


1.4  0 

0.63  -0.085 


Ki\ 


'-0.5  -0.5" 
-0.5  -0.5 


,Kp  i 


—0.5  -0.5" 
-0.5  -0.5 


h 

(36) 


Equilibrium  Point  2: 

U\e2  =  0.7,  U2e2  =  16  (deg),  Xle2  =  0.8826,  X2e2  = 


(34) 
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0.3678,  Te4  =  38.155  (N),  a4  =  0.6473,  and 


FIGURE  2.  JetCat  SPT5  turboshaft  engine  setup  on  the  test  stand 


0.6263,  Tel  =  181.9711  (N),  a2  =  1.0822,  and 


*2  = 


-2.83  -0.0008  p 
1.20  -2.10  5 


Kh 


-0.4  -0.4 
-0.4  -0.4 


KP2 


1.14  0 

0.78  -0.054 

—0.4  -0.4 


Cp  -I 

,  <^2  -  1  5 


(37) 


-0.4  -0.4 


-0.85  0.032 
0.32  -0.64 


1.1  0 

0.17  -0.011 


,  CP  =  I, 


Ki4  = 


-0.2  -0.2 
-0.2  -0.2 


,  Kp4  = 


-0.2  -0.2 
-0.2  -0.2 


(39) 


Equilibrium  Point  5  (Idle): 

=  0.145,  U2e5  =  16  (deg),  xle5  =  0.295,  x2e5 


U\e5  =  U-i'O,  U2e5  =  lO  (UegJ,  X\e3  = 

161,  Te5  =  7.317  (N),  a5  =  0.3361,  and 


0 


-0.38  -0.0008 
0.26  -0.34 


0.7  0 

0.1  -0.0024 


Ki5  = 


-0.1  -0.1 
-0.1  -0.1 


,Kp5  = 


-0.1  -0.1 
-0.1  -0.1 


cf  =  /, 

(40) 


FIGURE  3.  Ap  (a)  components  as  functions  of  parameter  a 


Equilibrium  Point  3  (Cruise): 

Mie3  =  0.4685,  u2e3  =  16  (deg),  xle3  =  0.7264,  x2e3  = 
0.5,  Te3  =  70.5125  (N),  a3  =  0.8818,  and 


-1.9  0.061 
0.45  -1.1 


Bp  = 


1.57 

0.3 


Kh 


-0.3  -0.3 
-0.3  -0.3 


,  Kp3  = 


0 

-0.023 


,  Cl  =/, 


-0.3  -0.3 
-0.3  -0.3 


(38) 


Equilibrium  Point  4: 

Uie4  =  0.3,  u2e 4  =  16  (deg),  Xle4  =  0.5327,  x2e4  = 


Other  controller  parameters  are  £c  =  1,  7]c  =  3.  The  ele¬ 
ments  of  Ap(a)  and  Bp(a)  matrices  have  been  shown  as  func¬ 
tions  of  scheduling  parameter  a  in  figures  3  and  4.  In  this  sim¬ 
ulation,  the  scheduling  parameter,  a,  is  defined  to  be  the  Eu¬ 
clidean  norm  of  the  gas  turbine  engine  spool  speeds,  which  are 
the  plant  outputs  and  capture  the  engine  nonlinearities.  Piece- 
wise  linear  interpolation  has  been  used  to  compute  matrices  in 
between  the  available  linearization  matrices  of  each  pair  of  ad¬ 
jacent  equilibrium  points.  The  equilibrium  values  of  the  plant 
states  and  control  inputs  are  shown  in  figure  5  as  functions  of 
scheduling  parameter  a.  Piecewise  linear  interpolation  has  been 
used  to  compute  equilibrium  values  in  between  each  pair  of  adja¬ 
cent  equilibrium  points.  The  equilibrium  manifold  in  a  3D  space 
of  two  spool  speeds  and  fuel  flow  control  input  is  shown  in  fig¬ 
ure  6.  The  elements  of  control  matrices  Kp(a)  and  KL(a)  have 
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0 


0 


-0.2 
bf  -0.4 
-0.6 
-0.8 


-0.2 
c?  -0.4 
-0.6 
-0.8 


0 

-0.2 
J?  -0.4 
-0.6 
-0.8 


x 

k 

N, 

FIGURE  4.  Bp(a)  components  as  functions  of  parameter  a 


FIGURE  7.  Kp{a)  components  as  functions  of  parameter  a 


/ 

y 

s 

/ 

/ 

3  0.5  1 

1  1.5  2 

g  \ 

/ 

9 

y 

• 

/ 

y' 

3  0.5  1 

1  1.5  2 

16.5 
16 

15.5 
15 


0  0.5  1  1.5  2  0  0.5  1  1.5  2 


FIGURE  5.  Xe  (a)  and  ue(a)  as  functions  of  parameter  a 


o 

-0.2 

sf  -°-4 
-0.6 
-0.8 

C 

0 

-0.2 

sf  -°-4 
-0.6 
-0.8 


\ 
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x 
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/ 

X 

\ 

0 

-0.2 
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( 
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X, 

\ 

X 
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1  1.5  2 
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FIGURE  8.  Ki(a)  components  as  functions  of  parameter  a 


FIGURE  6.  Engine  equilibrium  manifold  in  3D  space  of  spool  speeds 
and  fuel  control  input 


been  shown  as  functions  of  scheduling  parameter  a  in  figures  7 
and  8.  Piecewise  linear  interpolation  has  been  used  to  interpo¬ 
late  Kp(a)  and  Ki(a)  using  the  predesigned  indexed  linear  con¬ 
trollers,  which  are  given  in  equations  (36)  to  (40). 

4.2  Simulation  Results 

To  show  the  stability  of  the  closed-loop  system,  40  differ¬ 
ent  (30  equilibrium,  and  10  non-equilibrium)  linearizations  have 
been  used,  to  solve  inequality  (35),  in  Matlab  with  the  aid  of 
YALMIP  [37]  and  SeDuMi  [38]  packages.  The  numerical  value 
for  the  common  matrix  P  is 

"  0.5232  0.0059  0.0913  -0.0177  -0.0293  -0.0011 ' 
0.0059  0.3406  0.0132  -0.0082  -0.0862  -0.0114 
0.0913  0.0132  0.1721  -0.0461  0.0044  0.0105 

-0.0177  -0.0082  -0.0461  0.1275  0.0388  0.0282  5 

-0.0293  -0.0862  0.0044  0.0388  0.2684  -0.0211 

-0.0011  -0.0114  0.0105  0.0282  -0.0211  0.2484 
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where  its  condition  number  is  6.8910.  Here,  we  implement  the 


time  (sec) 


FIGURE  9.  History  of  the  states  (. xp(t ))  for  the  nonlinear  system  and 
the  linear  parameter  dependent  model 


0  5  10  15  20  25  30  35  40 


time  (sec) 

FIGURE  10.  History  of  the  rate  of  states  (. xp(t ))  for  the  nonlinear 
system  and  the  linear  parameter  dependent  model 

proposed  parameter  dependent  gain  scheduled  controller  to  oper¬ 
ate  the  JetCat  SPT5  turboshaft  engine.  This  case  study,  simulates 
the  engine  acceleration  from  the  idle  thrust  to  the  cruise  condi¬ 
tion  and  then  its  deceleration  back  to  the  idle  condition  in  a  stable 
manner,  with  proper  tracking  performance,  for  the  standard  day 
sea  level  condition.  Simulation  results  are  shown  in  figures  9  to 
25.  Figure  9  and  10,  shows  the  history  of  the  nonlinear  system 
and  the  linear  parameter  dependent  model  states,  xp(t ),  and  the 
rate  of  states,  xp(t).  We  can  conclude  that,  the  linearized  model 
is  a  very  good  approximation  of  the  real  nonlinear  plant.  Fig¬ 
ure  1 1 ,  shows  the  history  of  the  norm  of  the  closed-loop  system 


FIGURE  11.  Norm  of  the  closed-loop  system  matrices  (||AC/||),  and 

(I  fell) 


matrices  |fe||,  and  |fe||.  The  figure  shows  the  boundedness 
of  these  two  matrices,  in  accordance  with  Assumption  1,  where 
k\  =  4.0327,  and  ks  =  2.1512.  Figure  12,  shows  the  history  of 
the  closed-loop  system  matrix  eigenvalues  A  {Aci }.  All  the  eigen¬ 
values  remain  negative  with  the  time  change  of  the  scheduling 
parameter  a.  Figure  13,  shows  the  history  of  the  schedul¬ 
ing  parameter  a  =  p(y)  =  \  \y\  \  =  \  \xp\  |  (the  Euclidean  norm  of 
the  engine  spool  speeds).  The  history  of  the  scheduling  param¬ 
eter  rate  a  =  x^px^ ,  also  has  been  plotted.  Both  a  and  a  are 
bounded.  Figure  14,  shows  the  phase  plot  for  both  spool  dynam¬ 
ics.  Figure  15,  shows  the  evolution  of  the  plant  states  which  are 
high  and  low  spool  speeds.  Figure  16,  shows  the  time  evolution 
of  the  controller  states. 

Figures  17  and  18,  show  the  outputs  (i.e.,  high  and  low 
spool  speeds)  tracking  their  reference  signals  properly.  Figure 
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FIGURE  13. 

(ry-  *PTXP  \ 

{<X-  IMF 


Scheduling  Parameter  (a  m  \  \xp  |  |)  and  its  rate  of  change 


FIGURE  15.  Plant  states:  high  and  low  spool  speeds  (. xp ) 


FIGURE  14.  High  and  low  spool  speeds  vs.  high  and  low  spool  ac¬ 
celerations 


19,  shows  the  history  of  the  thrust  and  it  is  following  its  refer¬ 
ence  command  from  idle  to  cruise  condition  and  then  back  to  the 
idle  for  standard  day,  sea  level  condition.  Figure  20,  shows  the 
evolution  of  the  control  inputs  v(t)  =  [vi(t),V2(t)]r ,  which  are 
inputs  to  the  augmented  system,  each  element  is  corresponding 
to  one  of  the  control  inputs  to  the  original  system.  Figure  21, 
shows  time  rates  of  fuel  and  prop  pitch  angle  inputs.  Figure  22, 
shows  fuel  flow  and  propeller  pitch  angle  histories  as  the  control 
inputs  to  the  plant.  Figures  23  and  24,  show  the  evolution  of  the 
controllers  integral  ( Ki(a ))  and  proportional  ( Kp(a ))  gain  matri¬ 
ces.  These  gains  have  been  obtained  by  interpolation  using  the 
predesigned  indexed  family  of  fixed-gain  controllers,  and  each 
controller  corresponds  to  one  equilibrium  point  of  the  engine. 
The  numerical  values  of  these  gains  are  given  in  equations  (38) 
to  (40),  which  represents  the  controller  gains  for  idle  and  cruise 
condition  and  one  more  equilibrium  point  in  between  these  two 


FIGURE  16.  Controller  states  (vc) 


FIGURE  1 7.  Output:  high  spool  speed  and  its  reference  signal 


operating  points.  Figure  25,  shows  the  histories  of  turbine  tem- 
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FIGURE  18.  Output:  low  spool  speed  and  its  reference  signal 


FIGURE  21.  Rate  of  change  for  fuel  and  prop  pitch  angle  control 
inputs 


FIGURE  19.  Thrust  and  its  reference  signal 


FIGURE  20.  Control  inputs  to  the  augmented  system  (v(f)) 


perature,  thrust  specific  fuel  consumption  (TSFC),  compressor 


FIGURE  22.  Fuel  and  prop  pitch  angle  control  inputs  (u(t)) 


pressure  ratio  and  corrected  air  flow  rate. 


4.3  Engine  Limit  Control 

To  handle  the  limits  on  the  turbine  engine  system  states  and 
control  inputs,  the  developed  gain  scheduled  control  system  can 
be  integrated  with  a  reference  governor.  Reference  governors 
have  been  developed  previously;  one  of  the  good  examples  of 
this  approach  is  presented  in  [41].  This  method  addresses  the 
problem  of  satisfying  input  and/or  state  hard  constraints  in  non¬ 
linear  control  systems.  The  approach  uses  receding  horizon  strat¬ 
egy  and  consists  of  adding  to  the  primal  compensated  nonlinear 
system  a  reference  governor.  The  proposed  reference  governor  is 
a  discrete-time  device  which  handles  the  reference  to  be  tracked 
in  an  on-line  fashion.  The  resulting  hybrid  system  satisfies  the 
constraints  as  well  as  stability  and  tracking  requirements  [41]. 
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time  (sec)  time  (sec) 


Comp.  Pressure  Ratio 


Corrected  Comp.  Airflow 


FIGURE  23.  Controllers  integral  gain  matrix  (Kj(a))  elements  histo¬ 
ries 


FIGURE  25.  Turbine  temperature,  TSFC,  compressor  overall  pres¬ 
sure  ratio  and  air  flow  rate  histories 


the  closed-loop  nonlinear  gas  turbine  engine  system  with  a  gain 
scheduling  controller.  Simulation  results  confirmed  the  applica¬ 
bility  of  the  proposed  controller  to  the  nonlinear  physics-based 
JetCat  SPT5  turboshaft  engine  model  for  large  transients. 
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FIGURE  24.  Controllers  proportional  gain  matrix  ( Kp(a ))  elements 
histories 

5  Conclusions 

A  MIMO  linear  parameter  dependent  model  of  the  nonlin¬ 
ear  gas  turbine  engine  system  has  been  developed;  and  a  gain 
scheduling  controller  with  stability  guarantees  for  this  system  has 
been  designed.  Piecewise  linear  interpolation  technique  has  been 
used  for  interpolating  the  parameter  varying  gain  scheduling  con¬ 
troller  in  between  the  predesigned  indexed  family  of  fixed-gain 
controllers.  The  scheduling  variable  in  the  design  process  is  an 
endogenous  parameter  (i.e.,  a  function  of  the  plant  outputs)  and  it 
has  been  defined  to  be  the  Euclidean  norm  of  the  gas  turbine  en¬ 
gine  spool  speeds.  Stability  of  the  closed-loop  gas  turbine  engine 
system  with  a  gain  scheduling  controller,  has  been  shown  using 
global  linearization  method.  It  also  has  been  shown  that  a  single 
quadratic  Lyapunov  function  can  be  computed  for  this  system  us¬ 
ing  convex  optimization  tools,  which  guarantees  the  stability  of 
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